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1.0 SUMMARY 


When 4D RNAV operations are used in a time-based en route metering environment, the 
performance of various descent strategies in mixed traffic operations and the compatibility of 
different descent strategies will become increasing areas of concern. An investigation was 
conducted of the effects on system throughput and fleet fuel of utilizing various descent 
strategies, both individually and in combination. The study was in support of the NASA 
Advanced Transport Operating Systems (ATOPS) Flow Management Avionics Research 
Studies program and was performed under NASA contract. The three strategies considered 
were clean-idle Mach/CAS, constant flightpath angle (CFPA) Mach/CAS, and fuel optimal. 

The flow management evaluation model (FMEM), used in a previous NASA contract to study 
the effects of various levels of time navigation (TNAV) fleet equipage on air traffic control 
(ATC) operations and fuel efficiency, was modified for this study. The model was updated to 
include path definition capabilities for three Boeing airplane types, the B737-300, B747-200, 
and B767-200, which were assumed to represent the gamut of air carrier turbojet traffic. 
Traffic inputs consisted of airplane pairs that entered the simulation at a common altitude 
(37,000 ft), common speed (0.78 Mach), and common distance (200 nmi) from the meter fix. 
Three 4D RNAV absorbable (elapsed) times were used. 

Results indicated that, for the assumed simulation conditions and arrival traffic distribution 
at a typical en route metering (ERM) airport where the 737 category of airplane comprised 
almost 88% of all arrivals, the optimal strategy had a fuel consumption advantage at the 
expense of a slight reduction in throughput over the other two strategies. However, parametric 
analyses indicated that the CFPA and, to a lesser extent, the clean-idle were less susceptible 
to throughput degradation than the optimal when airplane type mixes and strategy mixes 
were varied. 
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2.0 INTRODUCTION 


The NASA Advanced Transport Operating Systems (ATOPS) Flow Management Avionics 
Research Studies program is identifying requirements for the use of advanced avionics 
capabilities in the current and projected FAA time-based en route metering system (ref 1). 
General and interim path definition algorithms have been designed and an interim algorithm 
tested in simulation and in flight. 

The Federal Aviation Administration (FAA) has several advanced flow management concepts 
under study or development. These include the currently operational en route metering, the 
automated en route ATC, and integrated flow management concepts for the future. These 
programs and concepts are the basic elements of the National Airspace System Plan providing 
transition from today’s controller intensive ATC system to the highly automated system of 
2000 and beyond. 

The definition of alternative flow management concepts, the development of scenarios for 
ATOPS application, and the analysis and recommendation of further testing of scenarios were 
accomplished under contract NASl-14880 TRA-112. Additional work was performed under 
NASA Thsk Requirement 1-2 of contract NASl-16300. This included the application of a 
fast-time simulation of the Denver Air Route Traffic Control Center (ARTCC) to study the 
effects of various levels of TNAV equipage on ATC operations and fuel efficiency. The equipage 
study assumed all aircraft employed conventional clean-idle Mach/CAS descents. 

This report describes the results of an analysis of three descent strategies and their effects on 
throughput (number of aircraft processed per unit time) and fleet fuel. Three strategies were 
considered: (1) the clean drag and idle thrust descent using Mach/CAS speeds, (2) a constant 
flightpath angle descent with Mach/CAS speeds, and (3) the use of an optimal speed/thrust 
schedule to achieve ATC time control objectives. This study was defined under Task 6 of NASA 
contract NASl-17635. 

This report contains objectives, assumptions, and methodology employed, a description of the 
candidate descent strategies considered, and an analysis of the resultant traffic throughputs 
and fleet fuel consumed. Conclusions and recommendations are included. A mathematical 
description of the optimal strategy is included in Appendix A. 
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3.0 SYMBOLS AND ABBREVIATIONS 


Vsl 

■ns2 

I 

(j2 


% 

AAI 

ACPH 

ARTCC 

ATC 

ATOPS 

CAS 

CFPA 

Ad 

DOC 

^ef 

di 

d, 

ERM 

FAA 

FCLT 

ft 

Fk 

FMEM 

4D 

Ah 

hi 

ha 

JFK 

KCAS 

kn 

LAX 

NASA 

lb 

min 

nmi 

ORD 

Pcf 

Pci 

Popt 

P 737 

P 747 

P767 

Rm 

RNAV 


descent strategy probability of lead airplane 
descent strategy probability of trail airplane 
probability of descent strategy combination {m} 
analysis time interval 
throughput variance 
increase in time separation 
frequency of airplane combination{k} 
percent 

airplane arrival interval 
aircraft per hour 
air route traffic control center 
air traffic control 

advanced transport operating systems 
calibrated airspeed 
constant flightpath angle 
horizontal separation between aircraft 
direct operating cost 

minimum ATC horizontal separation, 5 nmi 

distance of lead airplane from meter fix 

distance of trail airplane from meter fix 

en route metering program 

Federal Aviation Administration 

freeze calculated landing time 

feet 

maximum flow rate of airplane combination {k} 
flow management evaluation model 
four dimensional 

vertical separation between aircraft 

true altitude of lead airplane 

true altitude of trail airplane 

John F. Kennedy International Airport 

calibrated airspeed in knots 

knot 

Los Angeles International Airport 

National Aeronautics and Space Administration 

pound 

minute 

nautical miles 

Chicago O’Hare International Airport 
frequency of CFPA strategy aircraft arrivals 
frequency of clean-idle strategy aircraft arrivals 
frequency of optimal strategy aircraft arrivals 
frequency of B737 airplane arrivals 
frequency of B747 airplane arrivals 
frequency of B767 airplane arrivals 

maximum flow rate for mixed strategy traffic combination {m} 
radio navigation 
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3.0 SYMBOLS AND ABBREVIATIONS (Concluded) 


sec second 

T system throughput 

At)5 minimum time separation for airplane combination {k} 

(Atk)i initial time separation for airplane combination {k} 
TNAV time navigation 
Vg2 groundspeed of trail airplane 
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4.0 OBJECTIVES, ASSUMPTIONS, AND METHODOLOGY 


The objectives of the study were to analyze and recommend a preferred descent strategy or 
strategies for use by 4D RNAV-equipped airplanes in future ATC operations that will have 
time-based metering. The strategy options evaluated for this study were clean-idle Mach/CAS 
descents, constant flightpath angle Mach/CAS descents, and point-mass optimal descents 
using variable speed and thrust schedules throughout the descent. The analysis and 
recommendation of one strategy over another were based on two criteria; fleet fuel usage and 
system throughput. 

The study assumed a future ATC environment when most air carrier and other high 
performance aircraft will have advanced flight management systems with time navigation 
(TNAV) capabilities. Also assumed for this future ATC environment was an advanced air 
traffic control system that will have the ability to maximize the traffic throughput rate by 
adjusting spacings over a meter fix based on predicted relative separation requirements. This 
ATC advanced flow management system would be an extension of the current en route 
metering (ERM) system. 

A composite U.S. fleet mix for the 1995 time period was assumed. This mix was based on 
extrapolations of current traffic trends at Denver, Dallas-Ft. Worth, and Minneapolis-St. Paul 
international airports. In the study, the effects of other fleet mix assumptions were also 
evaluated. Study results for Chicago O’Hare, John F. Kennedy, and Los Angeles international 
airports were also included. 

The approach that was taken in this study to determine a preferred descent strategy is 
illustrated in Figmre 1. The two steps taken for this analysis were first the performance of a 
separation analysis and then a throughput analysis. The separation analysis determined the 
minimum separation required for zero conflicts between a sequential airplane pair, flying 
from 200 nmi and cruise altitude to a meter fix at 10,000 ft and 250 KCAS. For each airplane 
in the traffic pair, inputs to the separation analysis included an assumed elapsed time, descent 
strategy, airplane t 3 rpe, and weight. Both airplanes entered the simulation at a common 
altitude of 37,000 ft and speed of 0.78 Mach in order to examine the influence only of descent 
strategy on conflicts (and therefore throughput), while controlling the effects of differences in 
initial altitude and speed. The pair was initially separated at cruise altitude by 5 nmi, which 
corresponds to a time separation of 0.67 min. Outputs of the separation analysis included 
zero-conflict fuel usage for each airplane and closest approach separation for the pair. For each 
airplane pair, a maximum flow rate value was computed. Maximum flow rate is the 
equivalent number of airplanes crossing the meter fix, while satisfying the zero-conflict 
criterion. These individual maximum flow rates for aircraft pairs were inputs to the 
throughput analysis. 

The throughput analysis generated an expected flow rate for a particular arrival traffic 
distribution by airplane type. Required inputs were traffic frequencies, i.e., the probabilities of 
occvurence in the surival fleet of each of the pairs. The throughput analysis combined the 
maximum flow rate and fuel data for each pair with the frequency data to determine a fleet 
fuel usage and the system throughput. The throughput analysis was performed for each of the 
three strategies over a parametric range of airplane type mixtures (at all three elapsed times) 
and over a parametric range of descent strategy mixtures for an assumed typical arrival 
distribution (for the short elapsed time only). 
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5.0 MODELING 


A revised version of the flow management evaluation model (FMEM), originally 
developed by Boeing to conduct a NASA-funded investigation of mixed TNAV equipage traffic 
in an en route metering (ERM) environment (ref. 2), was used to conduct the sensitivity 
analysis. 


5.1 FMEM ADAPTATION 

The FMEM simulates arrival traffic operations at an airport participating in the ERM 
program. The model was specifically configured to represent the Denver Air Route Traffic 
Control Center (ARTCC) airspace and arrival operations at Denver Stapleton International 
Airport. Airplanes are assigned landing slots based on a first-come first-served mechanization 
and the minimum time separation interval, known as the airplane arrival interval (AAD. The 
resultant schedule determines required meter fix times for all arrival traffic. ERM 
guarantees, or “freezes,” an airplane’s landing time when it is within a fixed number of flying 
minutes from its expected meter fix arrival time. This fixed number of minutes is called the 
freeze calculated landing time (FCLT) parameter. A TNAV-equipped airplane is presumed to be 
given its meter fix time assignment when it enters the freeze region at its freeze time. The 
difference between the airplane’s assigned meter fix time and its freeze time is its required 
elapsed time, which the airplane must absorb to make good its meter fix time. The assignment 
process is a dynamic one and can produce varying elapsed times, depending on traffic demand. 

The FMEM required adaptation to conduct an evaluation of descent strategies. Two 
additional descent strategies (CFPA and optimal) were added and a different algorithm for 
computing surveillance data was included. Some tabular performance data were replaced by 
polynomial expressions. Finally, the capabilities of computing traffic flow rate and inputting 
elapsed times were new requirements. 

5.1.1 Additional Descent Strategies 

Constant flightpath angle (CFPA) and fuel-optimal descent algorithms, both developed by 
Boeing for other studies, were adapted and included as two additional descent strategies for 
use by TNAV-equipped airplanes. These are described in greater detail in Section 6.0. 

5.1.2 Surveillance Calculations 

The optimal descent logic computes energy states in relatively small, discrete altitude steps 
(500 ft) in descent and time intervals (10 sec) in cruise. Linear predictions within these 
intervals will provide more accurate position and energy states than the corresponding 
waypoint-referenced calculations used in the CFPA and clean-idle implementations. The 
assumption of linear performance characteristics between wa 5 q)oints, which can be spaced 
considerably far apart, can lead to errors in predicting an airplane’s position and energy state. 
It is important to minimize these errors when path profile data of two nearby aircraft are 
compared to determine their relative spacing for conflict detection. While both the CFPA and 
clean-idle profile computation processes have always used these same integration step sizes, 
the intermediate calculation results were never saved but used only to compute 
waypoint-referenced data. The FMEM was adapted and these intermediate data are now the 
basis for calculating present position data. 
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5.1.3 Polynomial Expressions for Performance Data 

The optimal logic requires smooth polynomial expressions for drag polar and fuel flow. Since 
the optimal logic was initially developed for the B737-300 and B767-200, the performance data 
for these airplane types were already in polynomial form. Therefore, fuel flow and drag polar 
expressions were required only for the B747-200. All other performance data were not 
required to be in polynomial form, which include speed limit and thrust data. 

5.1.4 Elapsed Time Inputs 

As described above, the assignment of meter fix times in the baseline FMEM version is a 
dynamic process and depends on traffic demand. The current analysis required that all 
aircraft absorb the same elapsed time. The metering logic thus required adaptation to 
accommodate the assignment of ftxed elapsed times. 

5.1.5 Computation of Flow Rate 

The implementation of the flow rate calculation is discussed in the context of the deflnition of 
flow rate in Section 7.2. 


5.2 TRAFFIC MODEL 

The traffic inputs consisted of three airplane t 3 q>es. Each t 3 q)e utilized two representative 
gross weights, which were designated as “light” and “heavy.” Historical 1984 arrival traffic 
data at three ERM airports were consulted to create an average annual distribution by 
airplane type. These distributions were the bases for the weighting factors applied to the 
sensitivity analysis results to obtain typical throughput figures as a function of descent 
strategy and elapsed time. 

5.2.1 Airplane Types and Weights 

Assumptions were made to equate all other commercial turbojet airplane types to the 
B737-300, B767-200, and B747-200. The substitutions are summarized in Table 1 and were 
justified on the basis that older airplane types will most likely be replaced by their modern 
equivalents. These three t 3 q)es were also considered to nm the gamut of aircraft weights and 
performance. Two weight categories (“light” and “heavy”) have been assigned to each 
airplane t 3 q>e. The selection of the weight range was dictated by two considerations: (1) a 
realistic range of approach weights and (2) a parametric compromise between maximum 
weight range and maximum delay margin. The performance characteristics of the B737-300 
(CFM56-3-B1 engines), B767-200 (JT9D-7R4D engines), and B747-200 (RB-211B engines) were 
modeled. 

Each combination of airplane type and weight for either the clean-idle or CFPA descent 
strategy, and initial and terminal conditions, produced high-speed and low-speed (without path 
stretching) descent times, corresponding to speed schedules consistent with performance 
envelope limitations. Thble 2 lists the high- and low-speed (Mach/CAS) descent times for all 
combinations of airplane types, two weight categories, and two descent strategies. The optimal 
strategy was assumed not to produce the constraining descent times because that strategy, 
which is not limited to a Mach/CAS speed schedule, makes greater use of its performance 
envelope. This assumption was confirmed by subsequent optimal strategy runs. No airplane 
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was allowed to path stretch or fly faster than its high-speed descent time. Thus, certain 
airplane t 3 rpe and weight combinations deflned the fastest of all low-speed times and the 
slowest of all high-speed descent times. These comprised the short and long elapsed times, 
which were imposed on all airplane t 3 T>es and weight categories. These times necessarily 
meant that certain airplane t 3 T)es did not use their full-time envelopes. 

A higher B737 heavy weight would have produced a slower standard high-speed time. 
Similarly, a higher B747 heavy weight would have defined a faster standard low-speed time. 
Table 2 shows that the standard high-speed descent time of 1658 sec is determined by the 
heavy B737 (CFPA descent), and the low-speed time of 1819 sec by the heavy B747 (CFPA). 

5.2.2 Arrival Traffic Distribution 

The three airports that are currently involved in the ERM program are Denver Stapleton 
International, Dallas-Fort Worth International, and Minneapolis-St. Paul International. 1984 
annual arrival statistics were obtained (ref. 3) for these three airports. 

The U.S. domestic and international arrival distributions among the three types are shown in 
Thble 3. The table indicates that the proportion of the B737 type is essentially the same at 
these three airports, while the mixes between the B767 and B747 vary. In addition, the B737 
represents a substantial percentage (about 90%) of the total arrival traffic. Table 3 also 
summarizes air carrier turbojet arrival operations at all domestic U.S. airports and indicates 
that the ERM typical distribution is very similar to the U.S. average. The ERM average type 
mix was used to illustrate system throughput performance at a typical airport. 
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6.0 DESCENT STRATEGY CHARACTERISTICS 


The term descent strategy, as utilized in this study, referred to the definition of the 
algorithms and procedures for predicting the arrival flightpath (limited to the vertical 
flightpath in this study) of an airplane from the entry point at cruise altitude to the meter fix 
at 10,000 ft. The descent strategy determined the manner in which accelerations/decelerations 
were performed, the speeds flown in cruise and descent, and the use of thrust and drag. The 
descent strategy was constrained to traverse a fixed distance from entry point to meter fix 
during a fixed time period. The three descent strategies compared in this study were 
clean-idle, CFPA, and optimal. 


6.1 CLEAN-IDLE DESCENT 

The clean-idle strategy was based on descending at idle thrust with flaps, speed brakes, and 
gear retracted. The fixed elapsed time from entry point to meter fix was achieved using a 
modified binary search technique iteratively to determine the required speeds. A 
one-dimensional search was obtained by defining the cruise Mach to be equal to the descent 
Mach. Furthermore, the descent speed was defined in terms of a Mach/CAS family of speeds, 
where for every Mach there exists only one value of CAS. This was accomplished by 
partitioning the descent speed envelope using a critical altitude model. This model 
mathematically described the variation of critical altitude (the altitude where a constant 
Mach descent transitioned to a constant CAS descent) with true airspeed. Thus, for each value 
of Mach selected by the search algorithm, a cruise Mach and descent Mach/CAS. were fully 
defined. 

For each iteration of descent Mach, the clean-idle strategy computed an acceleration or 
deceleration at cruise altitude, a cruise segment, a descent segment, and another acceleration 
or deceleration segment at meter fix altitude. The cruise distance was adjusted so that the 
total distance for all segments equals the required distance from entry point to meter fix. 
Accelerations or decelerations were performed at maximum climb thrust or idle thrust, while 
in level flight in a clean configuration. A cruise segment was then computed at the descent 
Mach. 

Descent computations were based on point-mass, steady-state equations of motion that were 
solved for small altitude steps during descent. Rate of descent was computed as a function of 
thrust (set to idle thrust) and drag utilizing a high-speed drag polar formulation combined 
with lift as a function of weight. The descent was integrated over the entire altitude range by 
s ummin g time, distance, and fuel values at each altitude step. 

Figure 2 plots altitude and distance versus time for simulations of B737 and B747 clean-idle 
descents. The entry point altitude is 37,000 ft and is 200 nmi from the meter fix at 10,000 ft. 
The elapsed time from entry point to meter fix is 1658 sec (27.6 min). A common time 
reference is used for both airplanes and shows that the B737 is initially 5 nmi (0.67 min at 
Mach 0.78) ahead of the B747. Differing cruise/descent Machs for the two airplanes result in 
the B747 gradually reducing separation from the B737 to zero. At 14.5 min, the B747 reaches 
top of descent, while the B737 continues cruise until 17.6 min. The B737 remains above the 
B747 until 24.7 min where it crosses below the B747, at which point it has increased 
separation back to 5 nmi ahead of the B737. 
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6.2 CFPA DESCENT 


This strategy was based on descent at a constant inertial flightpath angle. Since winds were 
neglected in this study, the descent was also at a constant airmass flightpath angle. The 
algorithms and procedures were identical to those of the clean-idle strategy described in 
Section 6.1, except that the descent segment itself was computed differently. 

Descent calculations were based on point-mass, steady-state equations of motion solved for 
small altitude steps during descent. The rate of descent at each altitude step was defined 
geometrically from the descent speed and the flightpath angle (which was a fixed value). The 
thrust required was then computed as a function of rate of descent and drag. Drag was first 
computed based on a high-speed drag polar that assumed speed brakes were retracted. If the 
resulting thrust was less than the minimum available (idle thrust), then the drag was 
recomputed based on extending speed brakes such that the required flightpath angle was 
maintained using idle thrust. Time, distance, and fuel values were summed over the descent 
as in the clean-idle strategy. 

Figure 3 plots altitude and distance versus time for simulations of B737 and B747 CFPA 
descents. All conditions except descent strategy are identical to those of Figure 2. The B747 is 
shown to reduce separation gradually from the B737 to zero. Tbp of descent is reached at 14 
miTi for both airplanes. Differing speed schedules cause the B747 to fall behind and above the 
B737 diiring descent. 


6.3 OPTIMAL DESCENT 

This strategy utilized Pontryagin’s Minimum Principle to assert that a fuel-optimal airplane 
trajectory for a fixed terminal time was actually an optimal direct operating cost (DOC) 
trajectory for a certain cost of time. 

Algorithms for the optimal strategy were based on solving the general system of differential 
equations of motion commonly employed for aircraft trajectory computations in the 
longitudinal plane. Since the equations were fifth-order and nonlinear, a tenth-order, two-point 
boundary value problem resulted from application of Pontryagin’s Minimum Principle. 
Singular perturbation theory was utilized to simplify the solution by separating the system 
dynamics into slow and fast modes to reduce a high-order problem to a series of lower order 
problems. The resulting formulation separated the solution into the inner region (descent) and 
the outer region (cruise). Since the inner region still required solution of two-point boundary 
value problems, the energy state approximation was utilized to reduce this to only one state 
equation. A Fibonacci search parameter optimization scheme was then utilized to generate 
speed and thrust commands, which minimized DOC. DOC was defined in terms of time and 
fuel costs, which were combined to yield the cost index, defined as time cost divided by fuel 
cost. A Newton-Raphson search technique was then employed to find the cost index providing 
the required elapsed time. The resulting thrust, speed, and fuel parameters were optimal 
values. Further details of the mathematical formulation of this optimization are contained in 
Appendix A. 

Figure 4 plots altitude and distance versus time simulations of B737 and B747 optimal 
descents. All conditions except descent strategy are identical to those of Figures 2 and 3. The 
B747 is shown to reduce separation gradually from the B737 to zero at 7 min, and to begin 
descent at 10.5 min. The B737 begins descent at 16 min as the B747 moves ahead of the B737. 
The B747 then slows, and the B737 moves ahead at 23.3 min and drops below the B747 at 
24.5 min. 


10 



7.0 TRAFFIC FLOW RATE ANALYSIS 


This section addresses the question of what arrival flow rate can be sustained by each of 
the descent strategies under consideration. Assumptions were made that all aircraft were 
constrained to the same initial cruise altitude and speed, that the same 4D elapsed time 
requirements were imposed on all traffic, and that an advanced ATC system adjusted spacings 
at the entry point to ensure conflict-free operations throughout all descents. A separation 
analysis was conducted to determine minimum time separation data for all pairs of aircraft 
type and weight. Then, system throughput for each strategy and elapsed time was computed 
from the minimum time separation data and an assumed arrival traffic distribution. 
Additional throughput studies were made for three major U.S. airports that represented other 
arrival traffic distributions. Finally, the effect on system throughput of mixing descent 
strategies was examined. 


7.1 SEPARATION ANALYSIS 

Since the intent of the analysis was to determine minimum time separation (referenced at the 
simulation’s entry point) to support a minimum ATC distance separation throughout the 
descent, the computation of the loss in flow rate was done as part of a conflict detection 
analysis. That is, when a loss of minimum separation occurred, a determination was made of 
the corresponding loss in flow rate. An ATC separation violation was said to occur when two 
proximate airplanes were simultaneously within 5 nmi horizontally and 2000 ft (above FL290) 
or 1000 ft (at or below FL290) vertically of each other. The result of the separation analysis 
was a set of minimum time separations at the entry point for a given elapsed time and 
descent strategy. 

7.1.1 Simulation Assumptions 

Traffic inputs were pairs of aircraft, consisting of any two combinations of airplane type and 
weight category. Therefore, for each descent strategy and elapsed time, 36 runs were made. In 
addition to constraining each airplane to a fixed initial altitude and speed, each pair was 
initially separated by 5 nmi. This corresponded to an initial time separation of 0.67 min at 
37,000 ft and 0.78 Mach. Both aircraft also ended their descents at the meter fix at an altitude 
of 10,000 ft and speed of 250 kn (CAS). Since the aircraft arrival interval (AAI) was also set at 
this initial time separation and every airplane had to absorb the same elapsed time for a 
specific run, the pair was also separated by 0.67 min by the time the trail airplane reached the 
meter fix. Because of the differences in speed between the aircraft pair throughout the 
descent, time separation varied considerably. The analysis time interval | was defined as the 
interval between the freeze time of the first airplane and the meter fix arrival time of the 
second airplane. 

7.1.2 Separation Data 

Horizontal and vertical separation are defined as 

Ad - di - 
Ah — hi - hi 

where 1 and 2 refer to the lead and trail aircraft, respectively, d, is the distance (in nmi) of the 
ith airplane from the meter fix, and hj is the true altitude (in feet) of the ith airplane. These 
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data were computed at every simulation clock step throughout the analysis time interval 
Separation histories (Ad versus Ah) for the example pair of aircraft discussed in Section 6 are 
plotted in Figxires 5 through 7. Each plot presents separation histories for all three strategies 
and for a particular elapsed time, starting at an initial horizontal separation of 5 nmi (Ad = 5) 
and vertical separation of 0 feet (Ah = 0) and ending at a nominal Ad = 3.23 and Ah = 0. Each 
curve represents the progression of relative separation. Successive points on a curve are 
separated by a simulation time interval of 30 sec. 

In the figures, for both the optimal and clean-idle strategies, the B747 is shown to overtake 
the B737 (Ad < 0) at cruise altitude. Figure 5 is the descent speed profile plot for the short 
elapsed time case. In all three strategies, it is seen that initially the B737 decelerates at cruise 
altitude in order to meet its required elapsed time. Conversely, the B747 accelerates at cruise 
altitude. In the cases of the clean-idle and CFPA strategies, these characteristics are a result 
of the Mach/CAS speed schedules dictated by the respective descent strategy implementations. 
The behavior in the optimal strategy case arises when the two airplane t 3 rpes change to their 
respective speeds at cruise altitude to minimize fuel. The magnitudes of the changes for both 
airplane t3q)es are greater for the optimal strategy than those for the other two and explain 
the greater horizontal separation (at Ah = 0) excursions for the optimal in Figure 5. For the 
long elapsed time, the optimal aircraft exhibit the same characteristics, while in the CFPA 
and clean-idle cases, the B747 decelerates instead. 

The B747 begins its descent before the B737 (Ah < 0). The B737 eventually gains on the B747 
in altitude until it goes lower in altitude briefly before both aircraft reach the meter fix. In the 
CFPA case, the B737 loses its horizontal lead position to the B747 but descends first (Ah > 0), 
as Figure 5 shows. 

The magnitudes of the horizontal displacement in the regions where ATC separation criteria 
have been violated are indications of the relative increase in time separation required of the 
trail airplane at the entry point to preclude conflicts. 

If the B747 were the lead aircraft, the separation histories are almost the reverse of the 
situation above. For all three strategies, the B747 accelerates as soon as it begins its delay 
absorption (Ad > 5) and begins* its descent first. This is exemplified in Figure 8, which 
corresponds to a required elapsed time of 1819 sec. This figime also illustrates the fact that the 
final horizontal separation among all three strategies is not necessarily 3.23 nmi, because 
each profile computation of a particular airplane type and weight is accurate to within ± 1 sec 
of the required elapsed time. 

Examination of speed profile data of the heavy B737 and heavy B747 for the short elapsed 
time, for example, indicates the differences in true air speed throughout the descent. These 
are illustrated in Figures 9 through 11. For all strategies, the B737 immediately decelerates 
while the B747 accelerates after their respective freeze times (elapsed time = 0). The final 
cruise speeds for the clean-idle and CFPA cases are determined by the speed transition 
(Mach/CAS) model, while the cruise speed for the optimal strategy is selected primarily for 
best fuel consumption. 
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7.1.3 Minimum Time Separations 


The expression for required time separation is given by 

At, = (A^ + T (1) 

where At^ = minimum time separation between airplane pair {k} 

A(tk)i = initial time separation between airplane pair {k} 

T = increase in time separation between airplane pair {k} 


In these series of runs, (At^), = 0.67 min. The increase in time separation is the increment of 
time beyond the initial time separation by which the trail airplane is required to lag the lead 
airplane to ensure conflict-free operations. Therefore, t is correlated to a loss of horizontal 
separation when a conflict develops. If d^ is the minimum ATC horizontal separation 
criterion. Ad is the horizontal separation at an arbitrary time t, and v ^2 is the groundspeed (in 
knots) of the trail aircraft at time t, t (in minutes) is then given by 

T = 60 (d„f — Ad)/Vg2 

Because of the 5-mile rule. 


T = (300 - 60Ad)/Vg2 (2) 

It is seen from Equation (1) that the minimiun Atj^ to prevent conflicts throughout the analysis 
time interval $ occurs when t is maximum within |. This minimum At^ is referred to in this 
analysis as the minimum time separation. The program computes At^ whenever separation 
minima are violated and compares its value with a previously stored maximum value. If 
required, the new maximum will replace the previous value. This procedure is carried out 
until the trail airplane crosses the meter fix. Within the limitations that these calculations 
are made at 30-sec intervals, a value close to the true minimum time separation will have 
been computed. 


These values of minimum time separation will be used to determine maximum flow rates in 
Section 7.2. 


Ihble 4 contains the resultant minimum time separation data. 

7.2 FLOW RATE ANALYSIS 


As discussed above, this concept of flow rate was based on having assured that ATC separation 
criteria were not violated. Flow rate is therefore related to minimum time separation, as 
defined in Paragraph 7.1.3. 

7.2.1 Flow Rate 


If minimum time separation At|^ (in minutes) is the conflict-free time separation between a 
pair of aircraft, defined as pair k, on the same route, then its flow rate (in aircraft per hour, 
ACPH) is 


F, = 60/ At, (3) 

If the loss of time separation t between a particular pair of aircraft does not reach its 
maximum until the trail airplane reaches the meter fix, then Ad in Equation (2) is determined 
by the time separation between the aircraft and the groundspeed of the trail airplane at the 
meter fix. 
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Maximum flow rates equivalent to the minimum time separation data of Table 4 are 
presented in Table 5. 

Because of the implementation of the 4D logic in the FMEM, path solutions were accurate to 
within ± 1 sec. Thus, if the lead airplane were 1 sec early and the trail airplane 1 sec late, the 
greatest time spacing at the meter fix was 0.70 min (0.67 min + 2 sec), with a final spacing of 
3.38 nmi. Because the groundspeed of the trail airplane was 288.85 kn (250 kn CAS at 10,000 
ft altitude), the limiting (worst case) flow rate could have been as high as 59.69 ACPH, based 
on calculations using Equations (1) and (3). The conclusion is that all flow rate values, 
assuming no other perturbations, cannot exceed a rate of 59.69. This value served as a check 
on the validity of the results of Thble 5. If the profile calculation had no time error, then the 
nominal limiting throughput was 57.77 ACPH. Note that, for the pair of airplanes consisting 
of the same type and weight, flow rate values were 57.77 in most cases. This is because the 
path solutions were exactly the same for both airplanes, the path solutions incurred the same 
error, and therefore the relative time spacing was the nominal airplane arrival interval, or 
0.67 min. The limiting flow rate calculation was based on multiple traffic on a common 
arrival route. Nevertheless, for the same reasons mentioned, the flow rate for any pair of 
aircraft, even when arriving on different routes converging at the meter fix, is still limited to 
59.69 ACPH. 

Adjusted flow rate, or the degradation of flow rate, as a function of time is plotted for each 
elapsed time in Figures 12 through 14. The final flow rate value at the end of | is referred to 
as the maximum flow rate for a given aircraft pair, strategy, and elapsed time. All three 
strategies are presented on the same plot for comparison. The cases are the same as discussed 
in Section 6 and Paragraph 7.1.2. The initial flow rate at the entry point is always 89.55 
because the initial time separation is 0.67 min. Elapsed time is referenced to the first 
airplane’s freeze time, or the beginning of These figures show that loss of minimum 
separation occurs immediately. Flow rate degrades steadily for all three strategies until 
somewhere in the middle of descent. The conflict-firee, maximum flow rates are achieved 
between 10 and 17 min, depending on strategy. 

7.2^ Parametric Traffic Analysis 

System throughput is here defined as the expected flow rate. This definition presumes a traffic 
distribution and large sample populations. Statistical analyses can be conducted in which the 
frequencies of airplane t}^es are used. 

Weighted averages and standard deviations can be computed. If weighting factors or arrival 
probabilities of pairs of aircraft are given by (if k represents a pair of particular aircraft 
t 3 T>e and weight mix), then throughput, T, is given by 

r = f 0,^* (4) 

where M is the number of combinatory pairs of airplane type and weight, or 36, and F)j is 
given by Equation (3). Its variance is computed by 

= S - 71 ' <5) 

*=i 


14 


For each descent strategy and elapsed time combination, families of throughput curves were 
generated from Equation (4) as functions of a particular airplane type percentage with a 
second type as a parameter. Varying airplane type mixes subject to 

P737 + PlAl + 1 

produces throughput sensitivity relationships, where P737, P747, and P7g7 are the arrival 
probabilities of the B737, B747, and B767, respectively. Figures 15 through 17 show the effect 
of varying B737 proportions with the B767 proportion as parameter cross-plotted with those 
with the B747 as parameter, for all three elapsed times and a clean-idle strategy. 
Corresponding plots for the CFPA and optimal strategies are presented in Figures 18 through 
20 and Figures 21 through 23, respectively. The clean-idle, CFPA and, to a lesser extent, the 
optimal results show that, for fixed B747 levels, throughputs are relatively constant for all 
ranges of B737 and B767 mixes. However, for increasing levels of the B747, throughput favors 
a higher percentage of the B767 than that of the B737. 

For the optimal strategy at fixed B747 levels, throughput performance depends more on the 
ratio between the B737 and B767 than do the clean-idle and CFPA strategies. For both the 
short and medium elapsed times and fixed levels of the B747 from 0% to 50%, minimum 
throughputs are generated when the B737 constitutes 50% of all arrivals. For the long elapsed 
time, 40% B737s are required for minimum throughputs for 0% to 50% fixed B747 levels. This 
implies that for the optimal strategy, system throughputs for 50% B737 levels will be 
minimum regardless of the relative mixes between the B747 and B767. However, these 
reductions in throughput are not as significant as when the percentage of B747s is increased 
from zero to some intermediate level, for either fixed levels of the B737 or the B767. 

Figures 15 through 23 also demonstrate that increasing the levels of the B747 has more of a 
deleterious effect on throughput for the optimal strategy than for the CFPA strategy. Ihble 6 
summarizes the deterioration of throughput under two scenarios. The first considers what 
happens when, with the traffic always maintained at equal B737 and B767 levels, the 
numbers of B747s is increased from 0% to 40% of all arrivals; and the second, when the 
percentage of B747s is increased from 0% to 50% when no B767s are part of the traffic. 

For all strategies, increasing the 4D times generally decreases throughput for the same traffic 
mix. This effect is less pronounced for the CFPA strategy than for either the clean-idle or 
optimal strategy. For example, with no B767s and when the B737 and B747 are equally 
distributed, throughput degrades by 3.58%, 1.54%, and 2.72% for the clean-idle, CFPA, and 
optimal strategies, respectively, when the 4D time is increased from 1658 to 1819 sec. 
Similarly, when the B737/B747 ratio is 30/70, throughput reductions are 4.14%, 2.22% 
and 6.02%. 

While the best possible throughputs are the same regardless of strategy, the worst 
throughputs are achieved by the optimal strategy. Thble 7 summarizes the best achievable 
results and the conditions under which they occmr. The worst throughputs are shown in Ihble 
8. In the worst possible scenario, when all aircraft engage the optimal strategy and there are 
no B767s, increasing the percentage of the B747 from 0% to 70% of all arrivals (and therefore 
reducing the B737 percentage from 100% to 30%) would reduce system throughput by 15.1%. 

7.2.3 Typical Arrival Traffic Analysis 

The t3T)ical traffic distribution (par. 5.2.1) is influenced by the predominance of the B737 type. 
Weighting factors for the typical traffic distribution are shown in Table 9. 
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Weighted throughput means and standard deviations derived from Equations (4) and (5) are 
given in Table 10. Weighted means for fuel are also provided. 

As shown in Table 3, the B737/B737 combinations contribute 77.2% (19.3% x 4) of all aircraft 
pairs in typical arrival traffic. The average of the individual flow rates for the B737, extracted 
from Table 5, in their four weight combinations, and the average minimum time separations 
related to these flow rates, are summarized in Table 11. Based on the fact that the B737’s 
optimum cruise altitude and speed are essentially 37,000 ft and 0.78 Mach, respectively, these 
results indicate that, for almost 80% of all traffic combinations, throughputs will be the same, 
regardless of descent strategy and required elapsed time. Furthermore, an entry point arrival 
interval of almost 1.04 min is sufficient to generate a throughput of approximately 58 ACPH, 
with no ATC intervention. A separation analysis of the B737 pairs shows that maximum flow 
rate is insensitive to descent strategy because an initial separation of almost 1.04 min 
(dictated by the constraints at the meter fix) is of such a magnitude that variations in 
separation behavior among the three strategies are not significant. 

7.2.4 Analyses of Other Distributions 

While the previous section examined a typical airport, some major U.S. airports experience 
somewhat different mixes of arrivals. The cases of John F. Kennedy (JFK) International, Los 
Angeles International (LAX), and Chicago O’Hare (ORD) were considered to analyze 
throughput performance as a function of different arrival traffic distributions. The summary 
of 1984 distributions by airplane type in Table 3 shows the B737 aircraft type only constitutes 
42.6% of all arrivals at JFK because of offloading of short-haul operations to nearby municipal 
airports. LAX and ORD show increasing frequencies of the B737. 

Weighted throughputs as functions of elapsed time for the t 3 q>ical ERM airport and three 
others are plotted in Figure 24. From the trends demonstrated in Figure 24, it is apparent 
that system throughput decreases as a single airplane type becomes less dominant in the total 
traffic mix. Furthermore, the optimal descent strategy is susceptible to greater throughput 
degradation than either of the other two strategies as traffic becomes more equally distributed 
among the three airplane types. 

7.2.5 Mixed Strategy Analysis 

Since it is probable that airlines will choose different flight management system 
implementations of descent strategy, the effect on throughput of mixing descent strategies was 
examined. For a given elapsed time, system throughput is the weighted average of all aircraft, 
weight and descent strategy combinations: 

T='2KR„ ( 6 ) 

m 

where combination m refers to an aircraft pair of specific type, weight, and strategy 
combinations and R„, is the maximum flow rate associated with m. The expression consists of 
324 terms since either of the pair can assume one of 18 airplane t 5 rpe-weight-strategy 
conditions. The frequency is also related to the weighting factor Wk of Equation (4): 

K = <^kVslVs2 

where «k is the arrival probability of an aircraft pair of t 3 rpe-weight combinations k, Tjgj is the 
descent strategy probability of aircraft 1 and ijg 2 is the strategy probability of aircraft 2. 
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A parametric study, for only the short elapsed time case, was also performed. In a fashion 
similar to the development in Paragraph 7.2.4, families of throughput curves were generated 
from Equation (6), which are functions of a particular descent strategy with a second strategy 
as a parameter, subject to: 


Popt Pci Pcf ^ 

where p^p^, p^j, and p^ are the arrival probabilities of optimal-, clean-idle-, and CFPA-strategy 
aircraft, respectively. The throughput sensitivity relationships are plotted in Figure 25, with 
the optimal strategy as the independent variable. These curves correspond to the assumed 
t5^ical traffic distribution and the short elapsed time. The trends expressed in the figture 
indicate that the optinial strategy has a greater influence on throughput than any of the other 
strategies. In particular, increasing the percentage of aircraft (for the assumed distribution) 
employing optimal descents degrades system throughput faster than corresponding increases 
of either the CFPA or clean-idle-equipped aircraft. On the other hand, an increase in CFPA 
aircraft as a percentage of total traffic improves throughput for either a fixed level of 
clean-idle or optimal aircraft. When none of the aircraft uses clean-idle, changing the fleet 
from 0% CFPA to full CFPA equipage increases throughput by 2.7%, whereas for traffic 
containing no CFPA aircraft, throughput improves by 2.6% when the traffic mix changes from 
0% clean-idle to 100%. Generally, for a fixed level of optimal aircraft, throughput is insensitive 
to variations in CFPA/clean-idle mixes, as Figure 25 shows, but slightly favors the CFPA as 
the optimal percentage increases. 

Figure 26 shows the results of a similar analysis performed for JFK International where the 
arrival traffic consists mostly of B737 and B747 aircraft in equal numbers. System throughput 
is shown to be insensitive to mixes between the clean-idle and CFPA strategies for a fixed 
level of optimal. Throughput appears to be a function primarily of the level of 
optimal-equipped aircraft. 


7.3 FUEL ANALYSIS 

In this analysis, the notion of maximum flow rate for an airplane pair assumed that the fuel 
used by the pair was independent of their relative spacing for a specific descent strategy and 
elapsed time. The requirement that no conflicts occur was satisfied by the model’s 
determining a sufficient minimum initial spacing rather than by issuing an ATC clearance. 
Therefore, the model computed for each airplane type the minimum achievable fuel 
consumption for a given descent strategy and elapsed time combination. If the minimum time 
separation (par. 7.1.3) between an airplane pair were reduced to increase system throughput, 
techniques would have to be employed to avert conflicts (speed reduction, nominal path offset, 
etc.). If any airplane were still required to make the same elapsed time, the loss in time 
related to path changes must be made up elsewhere in descent by a speed increase or, if 
possible, path shortening. These perturbations produce fuel penalties, which are not 
quantified here. As long as the system was able to support the required time separations 
without traffic demand exceeding throughput capacity, these zero-conflict fuel usage results 
could be considered applicable. 

Delay time was absorbed in each of the strategies by using a descent speed schedule which 
was unique for each elapsed time/strategy combination. The cruise portion of the descent 
accounted for 50% to 65% of the 200-nmi distance flown and was the major contributor to fuel 
usage. Comparisons of fuel usage were made for three values of delay times, 1658, 1739, and 
1819 sec. 
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7.3.1 Airplane Comparison 

Fuel bum values are shown in Table 12 for each airplane type, elapsed time, and strategy. The 
overall trend is for fuel usage to decline as elapsed time is increased, which is a result of the 
airplane being able to cruise at a lower, more fuel-efficient speed. 

A comparison of the differences in fuel usage can be seen in Table 13. Clean-idle is used as the 
baseline for each airplane and elapsed time. The optimal descent strategy reduced descent fuel 
burn by approximately 0% to 4%, while the CFPA strategy produced a penalty of as much 
as 6%. 

7.3.2 Fleet Comparison 

A combined weighted fuel burn has been derived to represent the typical mix of airplanes 
arriving at a typical ERM airport. The data, as in previous developments, were derived from 
the 1984 Airport Activity Statistics and used Dallas/Ft. Worth, Denver, and Minneapolis 
airival statistics. 

The fuel usage data are shown in Figure 27 for each strategy as a function of elapsed time. 
The data show that when the system throughput is under the limiting throughput value, the 
optimal strategy has the potential for a fuel advantage over the clean-idle strategy of 
approximately 1.5%, compared with a 2.6% penalty for the CFPA strategy. Fuel usage 
comparisons are also made for Chicago O’Hare, Los Angeles, and John F. Kennedy 
international airports (figs. 28 through 30) and demonstrate similar trends as for the 
typical airport. 
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8.0 CONCLUSIONS AND RECOMMENDATIONS 


The annual system throughput data in Table 10 assumed three descent strategies, three 
levels of absorbable (elapsed) times, specific initial traffic conditions, and a typical arrival 
aircraft type traffic distribution, in which the B737 comprises 88% of all arrivals. The CFPA 
descent strategy enjoys a slight throughput advantage over the clean-idle and optimal 
strategies. The differences among the strategies, however, are small for all elapsed times. This 
phenomenon results from the large numbers of the B737 airplane t 3 rpe in the arrival traffic as 
a whole. System throughputs, in fact, are close to the limiting throughput of 57.77 ACPH, as 
derived in Paragraph 7.2.1. The implications are that the limit will only be achieved if all 
traffic consisted of the same airplane type and weight combination, and that traffic composed 
of a more equitable distribution of airplane types would only decrease the expected 
throughput. If the traffic distribution were to show a smaller contribution by the B737 t 3 q>e 
and correspondingly larger influence of the other two types, system throughput rates for all 
strategies would deteriorate, with the optimal strategy incurring the greatest loss and the 
CFPA the least. These trends are demonstrated by the results of the Los Angeles, Chicago 
O’Hare, and John F. Kennedy international airports’ studies and the throughput sensitivity 
relationships in Figures 15 through 23. However, in the near- and far-term traffic 
environment, the actual arrival mix is forecast to remain similar to the t 5 q)ical distribution 
assumed for the analysis. 

Fleet-fuel savings results corresponding to conflict-free operations show the reverse trend. 
While the optimal strategy had the lowest throughput, it also used the least fuel. If the 
optimal strategy is used as the reference, fuel bum penalties of the clean-idle and CFPA 
strategies are 37.7 and 104.7 lb for the short; 17.1 and 58.7 lb for the medium; and 27.4 and 
57.6 lb for the long elapsed times, respectively. 

Therefore, the assumptions and conditions employed in this analysis favor the optimal 
strategy for fuel usage, with a slight disadvantage in throughput relative to the other two 
strategies for the assumed traffic distribution. However, if the individual maximum flow rate 
figures of Thble 5 bear some relationship to, or are in some way measures of, the extent to 
which ATC must intervene in order to increase throughput, then a greater fuel penalty will be 
experienced by those airplane combinations having lower maximum flow rates. Actual fuel 
penalties due to reduced flow rate cannot be established until conflicts can be dynamically 
resolved in the simulation. However, an approach is presented in Appendix B to quantify the 
fuel usage and conflict workload sensitivities to various arrival rates under the same distance 
and conflict-free assumptions made earlier. 

The preceding analysis assumes restricted airspace capacity because of the single (common) 
route and zero-conflict criteria. No excess delay absorption techniques or conflict resolution 
were considered within the system. Airspace capacity can be increased if traffic can be 
distributed over several paths. These distributed configurations include parallel paths on the 
same route, additional routes, or combinations of merging routes. These formulations enable 
the input (entry point) flow rates to match or exceed the output (meter fix) flow rate. However, 
because of separation requirements at the meter fix, its airspace feeder system throughput 
would still be constrained to the nominal limiting throughput value. Since system throughput 
with a typical ERM arrival traffic distribution is close to the limit, it would appear that an 
effort to analyze the effects of increasing airspace capacity would not justify the slight possible 
gains in throughput. As long as one airplane type dominates the arrival statistics, system 
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performance may be largely insensitive to descent strategy. When system throughput is 
expected to be near the limit, the study results suggest that no ATC or airline benefits will be 
achieved by increasing airspace capacity or decreasing the minimum, zero-conflict spacing 
required between aircraft. 

Implicit in this zero-conflict, single-route analysis is the assumption that the initial spacing 
between consecutive aircraft be made as large as required to satisfy the zero-conflict criterion. 
In the present analysis, that spacing was only a little more than 1 min on the average, or an 
in-trail spacing of 7.74 nmi. This also suggests that descent strategy is not an issue when the 
demand on the system is low. 

lb summarize, the comparative findings of the current analysis are as follows: 

1. For a typical arrival distribution, the optimal strategy is the most fuel efficient with a 
slight throughput disadvantage. 

2. The CFPA strategy generates the greatest throughput for zero conflicts, regardless of 
arrival traffic distribution, while the optimal generates the least. 

3. The optimal strategy consumes the least fuel per aircraft for zero conflicts, regardless of 
arrival traffic distribution, while the CFPA consumes the most. 

4. In mixed strategy traffic, throughput is the most sensitive to the percentage of 
optimal-strategy airplanes. 

5. The clean-idle and CFPA strategies can be mixed interchangeably with negligible 
throughput differences. 

Two additional analyses are suggested for firrther investigation of descent strategies for use in 
a metered environment. One study assumes arrival traffic in which the three airplane types 
enter the metering system at their optimum altitudes and speeds. The second examines the 
effect of changing the metering freeze radius (freeze calculated landing time, FCLT) on 
throughput performance. 

Whether the current results can be extrapolated to throughput performance of traffic in which 
airplane types enter the metering system at their optimum speeds and altitudes is not 
apparent. It appears that the combinations of slower airplanes entering the system ahead of 
faster ones would produce improved maximum flow rates in this configuration, because the 
initial conflicts at cruise altitude noted in the current analysis would probably be eliminated. 
This kind of analysis also lends itself to a statistical evaluation like the one in Section 7, as 
opposed to requiring a Monte Carlo simulation. Initial separations can also be defined in this 
context that will result in zero conflicts. Similar parametric analyses can be carried out to 
study throughput sensitivities to mixed traffic and mixed strategies. 

The effect of varying the FCLT parameter also would be a significant study since the relative 
proportions of descent and cruise segments will alter maximum flow rate performance. 
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APPENDIX A. THE OPTIMAL DESCENT STRATEGY FORMULATION* 


Aircraft Equations of Motion 

The aircraft model used for the optimal strategy formulation assumed a point mass 
approximation, that is, state variables describing the vehicle attitude were either omitted or 
used as control variables. The energy-state approximation of the point-mass longitudinal 
model of the aircraft, in the presence of wind, was taken to be in conventional state variables 
as follows: 


tl 
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where x was the range, V the airspeed, Vy, the windspeed, m the mass, f the fuel flow rate, E 
the energy height, y the flightpath angle, T the thrust, and D the drag. The second term on 
the right hand side of Equation (A.3) is normally small and can usually be neglected, t is a 
small “singular perturbation” parameter that arises as a consequence of the particular 
aircraft dynamics and an appropriate choice of scaling the equations of motion. The energy 
height E was related to the altitude and the groundspeed through 

E = h + (V + Vj/2g (A A) 

The airspeed V and the thrust T are the control variables, varying within the limits: 

< r < (^.5) 

< K < (^.6) 

Both Vnjin and were functions of altitude and represented the controllability and 

structural and performance limitations on the aircraft. The airplane model also involved the 
fuel flow rate f(h, M, T), the drag polar Cd(Cl,M), which is modeled using multiple regression 
methods, and minimum and maximum thrusts T ,„i„ (M, h), T „,„„ (M, h). M denotes the Mach 
number, the drag coefficient, and Cl the lift coefficient. 

Performance Index 

The optimization problem was to steer the system. Equations (A.l) through (A.3), from an 
initial state (x;, mj, Ej) at tj to a final state (Xf, mf, Ef) at fixed final time tf so that the fuel spent 
is minimized. Equivalently, the expression 

7= \^Cjfdt {A.l) 

was minimized where Cf was the cost of fuel. 


* Adapted from Four^Dimensional Fuel-Optimal Guidance in the Presence of Winds, A. Chakra varty, 
Journal of Guidance, Control, and Dynamics, VoL 8, No. 1, January-February 1985 
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Pontryagin’s Minimum Principle 

The Hamiltonian for Equations (A.l) through (A.3) and (A.7) was 


H, = Cff + KiV + KJ - \J- 


+ Xf 


(F + VJiT - D) 
mg 


M.8) 


where X^, X„, and Xg were the range, mass, and energy adjoint variables respectively. 
Pontryagin’s minimum principle states that the Hamiltonian is minimum along an optimal 
trajectory. Furthermore, since the final time is fixed, and Hj is not an explicit function of time. 
Hi is constant along the optimal trajectory and given by 


min 

T,V 


{HA = K 


{A. 9) 


K has the units of cost per unit time, and if Ct = -K is selected. Equation (A.9) may be 
rewritten 


min 

T,V 


c, + (C^ - KV + K(V + VJ 

(V+ VJ(T- D) 
mg 


+ Xe 


= 0 


(> 1 . 10 ) 


along the optimal trajectory. It therefore reduces to a direct operating cost (DOC) optimization 
problem with free terminal time and cost parameters Cf and Cj. The balance between fuel and 
time costs was expressed by the so-called cost index 



(A. 11) 


The traditional units used in the cost index are Ct in dollars per hour and Cf in cents per 
pound. The 4D optimization problem was now reduced to finding the correct Ct or Cl, for 
assigned arrival time tf. If the nominal trajectory corresponded to maximum range cruise 
(MRC), for which Ct = 0, a fuel-optimal trajectory for any delay would correspond to a negative 
cost of time, contrary to usual DOC choices. 


Cruise Cost Function 


As e -► 0, the outer solution (according to singular perturbation theory) was reduced to 


mm 

h,V C, + (Cf - \Jf + K(V + VJ 
(T=D) L 

Using Pontryagin’s minimum principle, we further get 


= 0 


_ , _ mm 

~ h,V 


c, -H (Cf - K)f 
iv -h FJ 


T=D 


where the ratio to be minimized is the cruise cost function. 


(^. 12 ) 


(>1.13) 
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Climb/Descent Cost Function 


During climb/descent, the independent variable was defined as t = t/e for climb and a = (t - 
tf)/e for descent such that as e 0, t, and a were finite over climb and descent, respectively. 
Equation (A. 10) was therefore transformed into 


min 

T,V 


c, + (C^ - K,V + KiV + y.) 


+ X, 


iV+ VJ(T - D) 
mg 


= 0 


(A. 14) 


where Xn,c and corresponded to the values during cruise, lb minimize Equation (A. 14), the 
energy adjoint 


Xf 


min 

V,D<T^T„^ 

max 

F, T„,„^T<D 


C, + (C, - X„,)/ + X,,(K + FJ 
(7 - D) (F + FJ/mg 


E 


(.A. 15) 


The minimization was done to get the climb solution, and the maximization to get the descent 
solution. The ratio in Equation (A. 15), to be optimized at current energy E, was called the 
climb/descent cost function. 
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APPENDIX B. DESCENT STRATEGY PERFORMANCE SENSITIVITY 

TO ARRIVAL RATE 


The main body of the report described an analysis whose results were based on separation of 
en route 4D arrival traffic to maximize throughput, without requiring ATC conflict 
intervention. Throughput and fuel usage performance did not assume accumulation of en 
route or terminal area delays, but instead depended on the existence of an advanced metering 
system that could space arrivals at the entry point to satisfy the zero-conflict criterion. 

On the other hand, if the rate of arrival traffic were increased, delays required to maintain 
minimum separation could accumulate to an extent that fuel usage figures for traffic 
conditions (airplane type mix and descent strategy) demonstrating lesser zero-conflict 
throughput performance will produce worse overall fuel economy. The analysis below attempts 
to quantify throughput, average fuel usage, and conflict workload for a range of input arrival 
rates for equally-spaced traffic. 

The basic assumption is that, if Atj is the minimum time separation between airplane pair i 
(interarrival time required to maintain nonconflicting separation between pair i; see par. 
7.1.3), the trail airplane will require delay of magnitude (At; - T) when input traffic time 
separation, T, is closer than Atj. That is, at some point, input arrival rate, AR = 1/T, exceeds 
the ability of the system to sustain that rate at the meter fix. Another assumption is that 
excess delay will be taken by the trail airplane by flying a vector at its original cruise altitude 
and speed, beginning immediately after receiving its 4D clearance. 


If input traffic of arbitrary sequence of airplane types consists of n airplane pairs, then total 
delay, D„, can be written as 

D„ = d^iT) -h [d,{T) + ^2(7) + • • • 

+ d,{T) -H • • • + d„{T)] (R.l) 


where D„ = total delay for n airplane pairs 

T = arrival spacing between airplane pairs 


d,{T) 



when Atj < T 
when Atj > T 


(B.2) 


Note that term (Atj - T) is the delay that the trail airplane of pair i must take when initial 
spacing (input interarrival time) is insufficient to maintain minimum time separation needed 
to relieve conflict between the pair, and that for any given delay required of pair i, the same 
delay has to be taken by all subsequent pairs through pair n to satisfy the same conflict-free 
criterion. This cascading phenomenon penalizes aircraft appearing later in the traffic list. 


Ibtal fuel usage for n airplane pairs, F^, is given by the expression 
F„ = FBj + [FB 2 + ff2di(T)] -I- ■ • • -I- [FB„^.i 

+ [d^iT) +•• - + d„{T)]] (B.3) 

where FBjj = conflict-free fuel used by the kth airplane in meeting its 4D time 
ffk = fuel flow of kth airplane at cruise altitude and speed 


and dj(T) is given by (B.2). 


25 


Tbtal fuel is shown to depend on total accumulated delays of all n airplane pairs. 


Finally, average fuel usage (per airplane) is given by: 


F„ = 


n +1 

For n airplane pairs average meter fix throughput, R, is given by 


(BA) 


(B.5) 


where t^i = elapsed time between the meter fix crossings of the first and last airplanes in the 
traffic. 


In turn. 


where 


tel = gi(T) + giiT) + • • • + g„(T) 



when Atj s T 
when Atj > T 


(5.6) 


(5.7) 


Expression (B.5), evaluated with (B.6) and (B.7), implies that average throughput is l/T when 
the arrival rate is small enough that no excess delay is caused by any airplane pair. Also, 
throughput can no longer match arrival rate when delays are needed to avoid conflicts. In 
particular, tgi decreases hsrperbolically with increasing arrival rate until the first applicable 
airplane pair comes into conflict. For sxifficiently large arrival rate, throughput no longer 
depends on arrival rate, but only on minimum time separations. Therefore, while throughput 
reaches saturation (remains constant), system delay increases according to (B.l). 


Note that throughput expression (B.5) is an alternative definition to that represented by 
Equation (4) in Paragraph (7.2.2). 

A measure of controller workload is the number of clearances that have to be issued to 
maintain minimum traffic separation. This is simply the number of trail airplanes among n 
pairs that require delay (vector) for a given arrival rate. Conflict count, is therefore: 

Neon = c,{T) + C2(70 + • • • + c„{T) (5.8) 

where 

when A/, s T 
when Atj > T 

The form of expressions like (B.l), (B.3), (B.6), and (B.8) implies sensitivity to the sequence of 
airplane types in the traffic for a fixed arrival rate. They also depend on the number of 
airplanes in the traffic list, that is, on the number n. 

Significant performance characteristics of throughput, fuel usage, and conflict workload can 
be approximated if simplifying assumptions are made. Specifically, if traffic were imagined to 
consist of a single airplane type whose performance characteristics of minimum time 
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separation. At, and fuel flow, f f , are “representative” (weightings) of the traffic mix at a given 
airport, then 

A/, = A /2 = • • • = A/„ = A/ 

//l = //2 = ' • • = ffn = ff 

The term At is the inverse of system throughput, such as those in Thble 10 for a typical ERM 
airport. The term fX can be considered a weighted (according to the assumed traffic mix) fuel 
flow of the three airplane types used for this study. Expressions (B.4), (B.5), and (B.8) reduce to: 

Approximate Fuel Usage 

< T 

> T {B.9) 


= 




FB 


when At 


FB + (At - T)ff when At 


where 


r = — 

AR 


and 


FB = 


1 

n + 1 


n+I 


1=1 


which is an expression for average imdelayed fuel usage (for known strategy and 4D time), 
some example values of which are plotted for large n in Figures 27 through 30 in the report. 


Approximate Meter Fix Throughput 


R = 



Approximate Conflict Count 


N = 

** ^con 


/ /I + 1 


when At < T 
when At > T 


when At < T 
when At > T 


(S.IO) 


(B.ll) 


The following relative conclusions can be drawn from expressions (B.9) through (B.ll), based 
on the observation that 


(^Oopi > (^Oc/mn > (^OcFPA (^-12) 

Equation (B.9) indicates that at sufficiently large arrival rates fuel usage performance of the 
three descent strategies changes in relation to each other because the optimal strategy has 
lower throughput (causes the airspace to reach saturation faster) than either the clean-idle or 
CFPA strategy. Stated in another way, zero-conflict throughput rates as defined in this report 
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can be viewed as measures of susceptibility to conflicts and therefore the need to resolve them 
by taking additional delays (using more fuel). Although the optimal strategy requires less fuel 
to fly than clean-idle or CFPA when arrival rate is low, that benefit is soon lost as conflicts 
begin to occur because the optimal strategy creates conflicts at lower arrival rates than the 
others (eq. (B.ID). 

Equation (B.IO) shows that among the three strategies the optimal causes meter fix 
throughput saturation at the lowest arrival rate. 

Breakpoints occur when input arrival rate equals system throughput for each strategy, which 
can be derived by equating the two parts of either expression (B.9) or (B.IO) to each other. 

A Monte Carlo simulation was also conducted to determine throughput, fuel usage, and 
conflict workload performance averaged over 1000 runs, using all three descent strategies at 
both the ERM and JFK international airports. Medium elapsed time (1739 sec) results of 
maximum flow rates and zero-conflict fuel usage numbers were used. Each of the 1000 runs 
per airport and strategy consisted of 120 airplane pairs whose sequence was randomly 
determined, but whose airplane type composition was constrained to mixtures enumerated in 
Table 3. Note that, for other values of n, results would be dependent on (but not proportional 
to) the number of airplane pairs analyzed. 

Means of fuel usage, meter fix throughput, and conflict workload are plotted in Figures B.l 
through B.6 and are evaluations of expressions (B.4), (B.5), and (B.8), respectively. Maximum 
values of the standard deviation to mean ratio are shown in Table B.l. 

These results demonstrate sensitivity to the distribution of minimum time separations {AtJ. 
The fuel and conflict figures illustrate two major phenomena: the effects of worst-case time 
separation behavior between the B737 type followed by the B747 t 3 rpe and the minimum 
attainable zero-conflict time separation. Tkble 5 in the main body of the report indicates at 
what arrival rates each strategy will experience separation violations among its pairings. 
B737/B747 combinations interact around 30 ACPH for the optimal strategy and around 40 
ACPH for clean-idle and CFPA. The effect of greater numbers of the B747 type is evident 
when JFK plots are compared to 'ERM plots at these arrival rates. In addition, the greater 
sensitivity of the optimal strategy to B767/B747 and B737/B767 pairs is also manifested for 
arrival rates around the high-30s through the mid-40s. 

The limiting flow rate in the mid-50s characterizes the minimum time separation of most 
airplane pairs and is graphically depicted in the conflict and fuel plots by the sharp rise in this 
region and by saturation in the throughput plots. 

Fuel and throughput performance relative to the CFPA strategy is plotted in Figures (B.7) 
through (B.IO). Clean-idle tracks CFPA performance closely, staying within 2.5% of CFPA’s 
fuel requirements and within about 1% of throughput. While requiring less fuel for arrival 
rates not exceeding 30 ACPH, the optimal strategy uses up to 14% (ERM distribution) and 
44% (JFK distribution) more fuel than the CFPA for greater arrival rates. Likewise, the 
optimal strategy’s throughput performance degrades by up to 5% (ERM) and 13% (JFK) for 
arrival rates in excess of 30 ACPH. 
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Figure 1. Descent Strategy Analysis Study Approach 
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Table 1. Airplane Equivalents 


AIRPLANE TYPE 

BOEING EQUIVALENT 

OC-9 


B727 


B737 

B737-300 

BAC-in 


BAC-146 


B707 


OC-B 


A-300 

B767-200 

A-310 


B757 


B767 


OC-10 


L-ion 

B747-200 

B747 
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Table 2. Elapsed Time Envelopes 


AIRPLANE 

TYPE 

WEIGHT 
(Ibx 1000) 

DESCENT 

STRATEGY 

MINIMUM 
ELAPSED 
TIME (sec) 

MAXIMUM 
ELAPSED 
TIME (sec) 

B737 

90 

Clean-Idle 

1609 

1883 

B737 

100 

Clean-Idle 

1648 

1887 

B737 

90 

CFPA 

1617 

1900 

B737 

100 

CfPA 

1658 

1900 

B747 

475* 

Clean-Idle 

1452 

2042 

B747 

564 

Clean-Idle 

1477 

1831 

B747 

475 

CFPA 

1455 

2029 

B747 

564 

CFPA 

1482 

1819 

B767 

215 

Clean-Idle 

1524 

1972 

B767 

270 

Clean-Idle 

1529 

2004 

B767 

215 

CFPA 

1523 

1955 

B767 

270 

CFPA 

1527 

1982 
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Table 3. 1984 Arrival Jaffic Distributions (Percent) 


(AIRPORT) 

B737-300 

TYPE 

B767-200 

TYPE 

B747-200 

TYPE 

John F. Kennedy 

42.6 

11.5 

45.8 

Los Angeles 

67.9 

6.2 

25.9 

Chicago O'Hare 

80.0 

7.9 

12.1 

Atlanta 

82.4 

11.0 

6.6 

Denver Stapleton 

86.2 

7.8 

6.0 

Oallas/Ft. Worth 

88.2 

4.1 

7.7 

Minneapolis/St. Paul 

90.6 

0.8 

8.6 

Typical ERH Airport 

87.9 

4.8 

7.3 

U.S. Average 

88.9 

4.6 

6.5 


Data derived from Airport Activity Statistics of Certificated Route Air Carriers, Office of Management 
Systems (Federal Aviation Administration) and Office of Aviation Information Management (Research 
and Special Programs Administration), December 1984. 
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Figure 2. Clean-Idle Descent Profiles, Elapsed Time = 1658 sec 


33 



DISTANCE (nmi) ALTITUDE (ft) 




Figure 3. CFPA Descent Profiles, Elapsed Time = 1658 sec 
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VERTICAL SEPARATION (ft) 

Figure 5. Comparative Separation Data— Heavy B737 Followed by Heavy B747, 
Elapsed Time = 1658 sec 
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HORIZONTAL SEPARATION (nmi) 



-Sooo -6o'o0 -4o'oO -2(loO o' 20o'o 

VERTICAL SEPARATION (ft) 

Figure 6. Comparative Separation Data— Heavy B737 Followed by Heavy B747, 
Elapsed Time = 1739 sec 
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’ -7000 -4500 -2000 500 

VERTICAL SEPARATION (ft) 

Figure 7. Comparative Separation Data — Heavy B737 Foiiowed by Heavy B747, 
Eiapsed Time =1819 sec 
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HORIZONTAL SEPARATION (nmi) 



Figure 8. Comparative Separation Data— Heavy B747 Followed by Heavy B737, 
Elapsed Time = 1819 sec 
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Figure 11. Descent Speed Profiles, CFPA Strategy, 1658 sec 
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Table 4. Minimum Time Separations (min) 


AIRPLANE 

PAIR 


1 658 sec 


1739 sec 



1819 sec 


Lead 

737H 

Tran 

737H 

Qpt 

1.039 

C-I 

1.039 

CFPA 

1.039 

Opt 

1T039 

C-I 

1.039 

CFPA 

1.039 

Opt 

1.039 

C-I 

1.039 

CFPA 

1.039 

737H 

737L 

1.054 

1.041 

1.050 

1.074 

1.053 

1.047 

1.048 

1.042 

1.041 

737H 

747H 

1.909 

1.347 

1.363 

1.965 

1.477 

1.448 

1.879 

1.592 

1.491 

737H 

747L 

1.977 

1.050 

1.075 

1.935 

1.107 

1.058 

1.760 

1.128 

1.037 

737H 

767H 

1.326 

1.049 

1.061 

1.304 

1.041 

1.041 

1.247 

1.045 

1.040 

737H 

767L 

1.124 

1.047 

1.051 

1.075 

1.045 

1.034 

1.037 

1.043 

1.034 

737L 

737H 

1.024 

1.036 

1.027 

1.003 

1.024 

1.030 

1.029 

1.035 

1.036 

737L 

737L 

1.039 

1.039 

1.039 

1.039 

1.039 

1.039 

1.039 

1.039 

1.039 

737L 

747H 

1.899 

1.366 

1.364 

1.976 

1.494 

1.447 

1.983 

1.615 

1.491 

737L 

747L 

1.965 

1 .072 

1.076 

2.002 

1.127 

1.057 

1.880 

1.152 

1-.037 

737L 

767H 

1.313 

1.047 

1.049 

1.377 

1.026 

1.032 

1.387 

1.041 

1.038 

737L 

767L 

1.109 

1.044 

1.040 

1.153 

1.030 

1.025 

1.146 

1.040 

1.032 

747H 

737H 

0.999 

1.036 

1.013 

0.971 

1.086 

1.070 

0.978 

1.090 

1.057 

747H 

737L 

1.014 

1.038 

1.025 

1.006 

1.101 

1.079 

0.987 

1.093 

1.059 

747H 

747H 

1.039 

1.039 

1.039 

1.041 

1.039 

1.039 

1.199 

1.046 

1.048 

747H 

747L 

1.064 

1.014 

1.045 

0.983 

1.088 

1.082 

1.123 

1,084 

1.054 

747H 

767H 

1.034 

1.046 

1.035 

1.007 

1.088 

1.073 

0.988 

1.096 

1.059 

747H 

767L 

1.025 

1.044 

1.026 

0.974 

1.092 

1.066 

0.977 

1.094 

1.052 

747L 

737H 

0.973 

1.060 

1.006 

1.027 

1.037 

1.027 

1.023 

1.044 

1.041 

747L 

737L 

0.988 

1.063 

1 .018 

1.062 

1.052 

1.036 

• 1.032 

1.048 

1.043 

747L 

747H 

1.013 

1.063 

1.032 

1.137 

1.138 

1.075 

1.223 

1.205 

1.041 

747 L 

747L 

1.039 

1.039 

1.039 

1.039 

1.039 

1.039 

1.147 

1.039 

1.039 

747L 

767H 

1.008 

1.070 

1.028 

1.063 

1.039 

1.030 

1.032 

1.050 

1.043 

747L 

767L 

1.000 

1.068 

1.019 

1.030 

1.043 

1.023 

1.022 

1.049 

1.037 

767H 

737H 

1.004 

1.028 

1 .016 

1.002 

1.036 

1.036 

1.030 

1.033 

1.037 

767H 

737L 

1.019 

1.031 

1.028 

1.037 

1.051 

1.045 

1.039 

1.036 

1.039 

767H 

747H 

1.474 

1.225 

1.230 

1.575 

1.339 

1.355 

1.596 

1.395 

1.425 

767H 

747L 

1.392 

1.007 

1.049 

1.362 

1.038 

1.047 

1.382 

1.027 

1.034 

767H 

767H 

1.039 

1.039 

1.039 

1.039 

1.039 

1.039 

1.039 

1.039 

1.039 

767H 

767L 

1.030 

1.036 

1.029 

1.005 

1.043 

1.032 

1.028 

1.037 

1.032 

767L 

737H 

1.012 

1.030 

1.026 

1.035 

1.032 

1.043 

1.040 

1.034 

1.043 

767L 

737L 

1.027 

1 .033 

1.038 

1.071 

1.047 

1.052 

1.049 

1.038 

1.045 

767L 

747H 

1.594 

1.303 

1.248 

1.682 

1.416 

1.364 

1.764 

1.463 

1.429 

767L 

747 L 

1.609 

1.009 

1.059 

1.603 

1.034 

1.054 

1.627 

1.029 

1.041 

767L 

767H 

1.047 

1.041 

1.048 

1.072 

1.034 

1.046 

1.173 

1.040 

1.045 

767L 

767L 

1.039 

1.039 

1.039 

1.039 

1.039 

1.039 

1.039 

1.039 

1.039 
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Table 5. Maximum Flow Rates (ACPH) 


AIRPLANE 

PAIR 


1658 sec 


1739 sec 


1819 sec 


Lead 

Trail 

OQt 

C-I 

CFPA 

Ogt 

C-I 

CFPA 

Qet 

C-I 

CFPA 

737H 

737H 

57.77 

57.77 

57.77 

57.77 

57.77 

57.77 

57.77 

57.77 

57.77 

737H 

737 L 

56.80 

57.73 

57.99 

57.74 

57.67 

57.72 

57.25 

58.47 

57.76 

737H 

747H 

31.54 

44.50 

44.40 

30,47 

40.26 

41.23 

32.09 

37.44 

40.08 

737H 

747L 

30.52 

56.46 

56.65 

30.58 

54.21 

56.93 

34.09 

53.21 

57.01 

737H 

767H 

45.63 

57.55 

57.81 

46.02 

57.67 

57.68 

48.13 

58.40 

57.69 

737H 

767L 

54.01 

57.58 

57.95 

55.82 

57.63 

57.12 

57.84 

58.41 

57.12 

737L 

737H 

58.77 

57.81 

57.55 

57.80 

57.87 

57.82 

58.30 

57.08 

57.78 

737L 

737L 

57.77 

57.77 

57.77 

57.77 

57.77 

57.77 

57.77 

57,77 

57.77 

737L 

747H 

31.72 

43.86 

44.14 

30.33 

39.66 

41.16 

30.39 

36.79 

40.06 

737L 

747 L 

30.71 

55.30 

56.19 

29.58 

52.99 

56.79 

31.91 

51.82 

57.01 

737L 

767H 

46.11 

57.58 

57.59 

43.57 

57.78 

57.74 

43.26 

57.70 

57.70 

737L 

767L 

54.81 

57.61 

57.73 

51.98 

57.73 

57.17 

52.34 

57.71 

57.12 

747H 

737H 

58.17 

57.50 

57.50 

58.34 

57.62 

56.94 

58.76 

56.46 

57.66 

747H 

737L 

57.18 

57.47 

57.72 

58.31 

57.52 

56.89 

58.23 

57.13 

57.66 

747H 

747H 

57.77 

57.77 

57.77 

57.59 

57.77 

57.77 

49.10 

57.40 

57.26 

747H 

747L 

57.31 

57.23 

57.44 

57.24 

57.60 

56.79 

50.59 

56.76 

56.91 

747H 

767H 

58.05 

57.29 

57.54 

58.24 

57.52 

56.85 

57.41 

57.06 

57.59 

747H 

767L 

57.78 

57.31- 

57.68 

58.37 

57.48 

56.30 

58.84 

57.08 

57.01 

747 L 

737H 

58.64 

58.05 • 

57.84 

58.89 

57.79 

57.93 

58.65 

57.47 

58.55 

747L 

737L 

57.64 

58.01 

58.06 

58.86 

57.69 

57.87 

58.12 

58.16 

58.55 

747 L 

747H 

58.24 

57.72 

58.11 

55.22 

51.93 

55.21 

50.77 

49.31 

51.41 

747L 

747L 

57.77 

57.77 

57.77 

57.77 

57.77 

57.77 

52.32 

57.77 

57.77 

747 L 

767H 

58.52 

57.83 

57.88 

58.78 

57.70 

57.84 

58.14 

58.09 

58.47 

747L 

767L 

58.25 

57.85 

58.02 

58.91 

57.65 

57.27 

58.72 

58.10 

57.88 

767H 

737H 

57.89 

57.99 

57.73 

57.87 

57.87 

57.86 

58.27 

57.16 

57.85 

767H 

737L 

56.91 

57.96 

57.95 

57.84 

57.76 

57.80 

57.75 

57.84 

57.84 

767H 

747H 

40.71 

48.83 

48.80 

38.32 

44.26 

44.01 

38.20 

42.47 

41.91 

767H 

747L 

43.21 

57.71 

57.66 

43.04 

57.84 

57.70 

43.40 

57.45 

57.08 

767H 

767H 

57.77 

57.77 

57.77 

57.77 

57.77 

57.77 

57.77 

57.77 

57.77 

767H 

767L 

57.51 

57.80 

57.91 

57.90 

57.72 

57,20 

58.34 

57.78 

57.19 

767L 

737H 

58.15 

57.97 

57.59 

57.75 

57.92 

58.44 

57.70 

57.14 

58.44 

767L 

737L 

57.17 

57.93 

57.81 

57.71 

57.81 

58.39 

57.18 

57.83 

58.43 

767 L 

747H 

37.63 

45.98 

48.24 

35.61 

41.82 

43.98 

34.61 

40.50 

42.00 

767 L 

747L 

37.34 

57.69 

57.52 

36.81 

57.89 

58.28 

36.87 

57.44 

57.66 

767L 

767H 

58.03 

57.74 

57.63 

57.56 

57.82 

58.35 

51.16 

57.76 

58.36 

767L 

767L 

57.77 

57.77 

57.77 

57.77 

57.77 

57.77 

57.77 

57.77 

57.77 





Tkble 6. Throughput Reduction Due to Increasing B747 Levels 


ELAPSED 

TIME 

(sec) 

STRATEGY 

B747 = 0% 
B737 = 50% 
B767 = 50% 

B747 = 40% 
B737 = 30% 
B767 = 30% 

REDUCTION 

(PERCENT) 


Clean-Idle 

57.77 

56.17 

-2.77 

1658 

CFPA 

57.77 

56.32 

-2.51 


Optimal 

55.80 

51.76 

-7.24 


Clean-Idle 

57.77 

55.30 

-4.28 

1739 

CFPA 

57.77 

55.64 

-3.69 


Optimal 

55.67 

51.46 

-7.56 


Clean-Idle 

57.77 

54.75 

-5.23 

1819 

CFPA 

57.77 

55.42 

-4.07 


Optimal 

55.54 

50.56 

-8.97 

ELAPSED 

TIME 

(sec) 

STRATEGY 

B747 = 0% 
B737 = 100% 
B767 = 0% 

B747 = 50% 
B737 = 50% 
B767 = 0% 

REDUCTION 

(PERCENT) 


Clean-Idle 

57.77 

55.80 

-3.41 

1658 

CFPA 

57.77 

55.92 

-3.20 


Optimal 

57.78 

51.15 

-11.47 


Clean-Idle 

57.77 

54.62 

-5.45 

1739 

CFPA 

57.77 

55.27 

-4.33 


Optimal 

57.77 

50.89 

-11.91 


Clean-Idle 

57.77 

53.80 

-6.87 

1819 

CFPA 

57.77 

55.06 

-4.69 


Optimal 

57.77 

49.76 

-13.87 








Table 7. Best Achievable Throughputs for Given Strategy and Elapsed Time 


ELAPSED 

STRATEGY 

THROUGHPUT 

CONDITIONS 

TIME 


(ACPH) 



(sec) 






Clean-Idle 

57.77 

B747-0*, 

any mix B737/B767 

1656 

CFPA 

57.77 

B747-0X, 

any mix B737/B767 


Optimal 

57.77 

B747-0X, 

all B737 or all B767 


Clean-Idle 

57.77 

B747«0X, 

any mix B737/B767 

1739 

CFPA 

57.77 

B747-0*, 

any mix B737/B767 


Optimal 

57.77 

B747-0X, 

B737=100<, B767-0X 


Clean-Idle 

57.77 

B747-0*, 

any mix B737/B767 

1819 

CFPA 

57.77 

B747-0X, 

any mix B737/B767 


Optimal 

57.77 

B747-0X, 

8737=100%, B767-0X 


l&to/e 8. Worst Achievable Throughputs for Given Strategy and Elapsed Time 


ELAPSED 

TIME 

(sec) 

STRATEGY THROUGHPUT REDUCTION 
(ACPH) FROM BEST 

(PERCENT) 

CONDITIONS 



Clean-Idle 

55.80 

-3.41 

B767-0X, 

B737=50X, 

B747-50X 

1658 

CFPA 

55.92 

-3.20 

B767-0X, 

B737=50X, 

B747-50X 


Optimal 

51.15 

-11.46 

B767-0X, 

B737=50X, 

B747-50X 


Clean-Idle 

54.56 

-5.56 

B767-0X, 

B737-40X, 

B747-60X 

1739 

CFPA 

55.27 

-4.33 

B767-0X, 

B737=50X, 

B747=50X 


Optimal 

50.89 

-11.91 

B767-0X, 

B737=50X, 

B747-50X 


Clean-Idle 

53.66 

-7.11 

B767-0X, 

B737-40X, 

B747-60X 

1819 

CFPA 

54.94 


B767-0X, 

B737-40X, 

B747-60X 


Optimal 

49.06 

-15.08 

8767-OX, 

B737=30X, 

B747=70X 
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FLOW RATE (ACPH) 




FLOW RATE (ACPH) 



FLOWRATE (ACPH) 



ELAPSED TIME (min) 


Figure 14. Flow Rate Degradation— Heavy B737 Followed by Heavy B747, 
Elapsed Time = 1819 sec 
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767 RATIO 



767 RATIO 



Figure 16. Mixed Traffic Throughput, Ciean-idie Strategy, Eiapsed Time = 1739 sec 
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THROUGHPUT (ACPH) THROUGHPUT (ACPH) 


767 RATIO 



Figure 17. Mixed Traffic Throughput, Clean-Idle Strategy, Elapsed Time = 1819 sec 


767 RATIO 



Figure 18. Mixed Traffic Throughput, CFPA Strategy, Elapsed Time = 1658 sec 
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THROUGHPUT (ACPH) THROUGHPUT (ACPH) 



Figure 21. Mixed Jaffic Throughput, Optimai Strategy, Eiapsed Time = 1658 sec 



Figure 22. Mixed Traffic Throughput, Optimai Strategy, Eiapsed Time = 1739 sec 
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THROUGHPUT (ACPH) 



Table 9. Probabilities of Airplane Pairs at a Typical ERM Airport 



AIRPLANE 

TYPE 

B737 


B767 

B747 

Airplane 

Type 

Weight 

Light 

Heavy 

Light 

Heavy 

Light 

Heavy 

B737 

Light 

0.19329 

0.19329 

0.01064 

0.01064 

0.01590 

0.01590 


Heavy 

0.19329 

0.19329 

0.01064 

0.01064 

0.01590 

0.01590 

B767 

Light 

0.01064 

0.01064 

0.00059 

0.00059 

0.00087 

0.00087 


Heavy 

0.01064 

0.01064 

0.00059 

0.00059 

0.00087 

0.00087 

B747 

Light 

0.01590 

0.01590 

0.00087 

0.00087 

0.00131 

0.00131 


Heavy 

0.01590 

0.01590 

0.00087 

0.00087 

0.00131 

0.00131 
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Table 10. Throughput and Fuel Use at a Typical ERM Airport 


ELAPSED 

TIME 

(sec) 

STRATEGY 

THROUGHPUT 

(ACPH) 

STANDARD 

DEVIATION, 

THROUGHPUT 

(ACPH) 

FUEL USE 
(lb PER 
AIRPLANE 


Clean-Idle 

57.26 

2.44 

1850.1 

1658 

CFPA 

57.28 

2.41 

1917.1 


Optimal 

55.69 

6.76 

1812.4 


Clean-Idle 

57.03 

3.25 

1750.2 

1739 

CFPA 

57.16 

2.96 

1791.8 


Optimal 

55.65 

7.01 

1733.1 


Clean-Idle 

56.88 

3.80 

1718.0 

1819 

CFPA 

57.17 

3.19 

1748.2 


Optimal 

55.76 

6.58 

1690.6 
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Table 11. Average Throughput and Minimum Time Separation Data for the B737-300 Type in 
Typical Traffic 


STRATEGY 

ELAPSED 

TIME 

(sec) 

AVERAGE 

THROUGHPUT 

(ACPH) 

AVERAGE 
MINIMUM TIME 
SEPARATION 
(min) 


1658 

57.78 

1.038 

Optimal 

1739 

57.81 

1.038 


1819 

57.77 

1.039 


1658 

57.77 

1.039 

Clean-Idle 

1739 

57.78 

1.038 


1819 

57.77 

1.039 


1658 

57.78 

1.038 

CFPA 

1739 

57.77 

1.039 


1819 

57.77 

1.039 
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Table 12. Total Fuel Usage (lb) 


AIRPLANE TYPE 

B737-300 

B767-200 

B747-200 

Delay 

(sec) 

Strategy 

Light 

Heavy 

Light 

Heavy 

Light 

Heavy 


Clean- 

Idle 

1459.4 

1541.4 

2517.5 

2738.3 

5320.6 

5843.0 

1658 

Constant 

Flightpath 

Angle 

1550.3 

1593.3 

2633.6 

2775.0 

5347.3 

5830.7 


Optimal 

1435.7 

1510.0 

2505.5 

2747.4 

5147.9 

5645.7 


Clean- 

Idle 

1365.2 

1438.1 

2390.2 

2634.7 

5178.7 

5776.8 

1739 

Constant 

Flightpath 

Angle 

1415.1 

1468.3 

2433.9 

2712.7 

5170.9 

5880.9 


Optimal 

1364.4 

1431.6 

2382.2 

2597.4 

5055.2 

5547.8 


Clean- 

Idle 

1326.9 

1408.1 

2314.9 

2584.5 

5134.0 

5845.1 

1819 

Constant 

Flightpath 

Angle 

1352.3 

>421.8 

2353.8 

2722.1 

5163.9 

6059.1 


Optima! 

1316.1 

1387.0 

2316.0 

2552.8 

5021.7 

5610.9 
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Table 13. Percent Change in Descent Fuel Usage Relative to Clean-Idle Descent 


AIRPLANE TYPE 

B737-300 

B767-200 

B747-200 

Delay 

(sec) 

Strategy 

Light 

Heavy 

Light 

Heavy 

Light 

Heavy 


Clean- 

Idle 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1658 

Constant 

Flightpath 

Angle 

6.23 

3.37 

4.61 

1.34 

0.50 

-0.21 


Optimal 

-1.62 

-2.04 

-0.48 

0.33 

-3.25 

-3.38 


Clean- 

Idle 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1739 

Constant 

Flightpath 

Angle 

3.66 

2.10 

1.83 

2.96 

-0.15 

1.80 


Optimal 

-0.06 ■ 

-0.45 

-0.33 

-1.42 

-2.38 

-3.96 


Clean- 

Idle 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1819 

Constant 

Flightpath 

Angle 

1.91 

0.97 

1.68 

5.32 

0.58 

3.66 


Optimal 

-0.81 

-1.50 

0.05 

-1.23 

-2.19 

-4.01 





THROUGHPUT (ACPH) 



CFPA— TYPICAL ERM 

CLEAN— TYPICAL ERM 
CFPA— CHICAGO 



CLEAN— CHICAGO 

CFPA— LOS ANGELES 
OPTIMAL— LOS ANGELES 
CFPA— JFK, NY 


CLEAN— LOS ANGELES 


OPTIMAL— CHICAGO 


CLEAN -JFK, NY 


OPTIMAL— LOS ANGELES 


OPTIMAL— JFK, NY 


ELAPSED TIME (sec) 


Figure 24. Throughput Performance at Various U.S. Airports 
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THROUGHPUT (ACPH) THROUGHPUT (ACPH) 



Figure 25. Mixed Strategy Throughput, Short Eiapsed Time, ERM Typicai Airport 



Figure 26. Mixed Strategy Throughput, Short Eiapsed Time, John F. Kennedy 
International Airport 
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FUEL (lb) 


Table B.1. Maximum Fuel Usage and Throughput Dispersion Ratios (Standard 
Deviation/Mean, a/p.) 




CFPA 

CLEAN-IDLE 

OPTIMAL 


Fuel 

0.012 

0.013 

0.032 

ERH 


Throughput 

0.002 

0.002 

0.006 


Fuel 

0.025 

0.025 

0.046 

JFK 


Throughput 

0.007 

0.007 

0.016 
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ARRIVAL RATE (ACPH) 


Figure B. 1. Fuel Usage (Per Airplane), Typical ERM Airport Distribution (Average of 1000 
Runs; 120 Airplane Pairs) 


64 







THROUGHPUT (ACPH) 



Figure B.3. Meter Fix Throughput, Typicai ERM Airport Distribution (Average of 1000 Runs; 
120 Airplane Pairs) 
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Figure BA. Meter Fix Throughput, JFK Intemationai Airport Distribution (Average of 1000 
Runs; 120 Airpiane Pairs) 
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ARRIVAL RATE (ACPH) 


Figure B.5. Conflict \Mirkload, Typical ERM Airport Distribution (Average of 1000 Runs; 120 
Airplane Pairs) 
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ARRIVAL RATE (ACPH) 

Figure B.6. Conflict Workload, JFK International Airport Distribution (Average of 1000 Runs; 
120 Airplane Pairs) 
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FUEL DIFFERENCE (PERCENT) 
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ARRIVAL RATE (ACPH) 


Figure B.7. Fuel Usage Penalty Relative to CFPA, Typical ERM Airport Distribution (Average 
of 1000 Runs; 120 Airplane Pairs) 
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Figure B.8. Fuel Usage Penalty Relative to CFPA, JFK International Airport Distribution 
(Average of 1000 Runs; 120 Airplane Pairs) 
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THROUGHPUT DIFFERENCE (PERCENT) 


0.0 



Figure B.9. Throughput Gain Reiative to CFPA, Typicai ERM Airport Distribution (Average of 
1000 Runs; 120 Airpiane Pairs) 
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